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Abstract 

Performing a fit to all publicly available data, we analyze the extent to which the latest 
results from the LHC and Tevatron constrain the couplings of the Higgs boson-like state 
at ~ 125 GeV. To this end we assume that only Standard Model (SM) particles appear in 
the Higgs decays, but tree-level Higgs couplings to the up-quarks, down-quarks and vector 
bosons, relative to the SM are free parameters. We also assume that the leptonic couplings 
relative to the SM are the same as for the down-quark, and a custodial symmetry for the 
V = W, Z couplings. In the simplest approach, the effective Higgs couplings to gluons 
and photons are computed in terms of the previous parameters. This approach is also 
applied to Two-Higgs-Doublet Models of Type I and Type II. However, we also explore 
the possibility that the net Higgs to gg and 77 couplings have extra loop contributions 
coming from Beyond-the-Standard Model physics. We find that the SM p-value ~ 0.5 is 
more than 2<r away from fits in which: a) there is some non-SM contribution to the 77 
coupling of the Higgs; or b) the sign of the top quark coupling to the Higgs is opposite 
that of the W coupling. In both these cases p-values ~ 0.9 can be achieved. Since option 
b) is difficult to realize in realistic models, it would seem that new physics contributions 
to the effective couplings of the Higgs are preferred. 
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1 Introduction 



The recent discovery u\ v2\ of a new particle with properties consistent with a Standard 
Model (SM) Higgs boson is clearly the most significant news from the Large Hadron Collider 
(LHC). This discovery was supported by evidence for a Higgs boson found by the CDF and 
DO collaborations at the Tevatron [3] and completes our picture of the SM. However, the SM 
leaves many fundamental questions open — perhaps the most pressing issue being that the SM 
does not explain the value of the electroweak scale, i.e. the Higgs mass, itself. Clearly, a prime 
goal after the discovery is to thoroughly test the SM nature for this Higgs-like signal. 

The SM makes precise predictions for the production cross sections of the Higgs boson 
H (via gluon-gluon fusion (ggF), vector boson fusion (VBF), associated production with an 
electroweak gauge boson V = W, Z (VH) and associated production with a tt pair (ttH)), and 
its decay branching fractions into various final states (77, ZZ&, WW^*\ bb, and tt) as a 
function of its unpredicted mass The observation of the Higgs boson at the LHC is based 

and WWW !J[~ 
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decay modes. The Higgs boson mass, 



primarily on the 77 Qg, ZZ& 
Mh, is quite precisely measured to be in the 125-126 GeV range using the high resolution 77 
and ZZ& final states (ll[[l2]Q The evidence for the Higgs boson at the Tevatron is based 
principally on the bb decay mode [3 
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the observed enhancements being consistent with 
a large range of possible Higgs masses. 

With the measurements in various channels, a comprehensive study of the properties of the 
Higgs-like state becomes possible and has the potential for revealing whether or not the Higgs 
sector is as simple as envisioned in the SM. In particular it is crucial to determine the Higgs 
couplings to gauge bosons and to fermions as defined by the Lagrangian 



C = 9 



CwmwW.W^ + Cz 



cos 9 



ir 
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(1) 



where the Cj are scaling factors for the couplings relative to their SM values, introduced to test 
possible deviations in the data from SM expectations. In principle all the Cj are independent, 
in particular the Cp can be different for up- and down- type quarks and/or leptons. A significant 
deviation of any Cj from unity would imply new physics beyond the SM. 

While fits to various combinations of C/'s are performed by the experimental collaborations 
themselves 



12,15 



we find it important to develop our own scheme in order to bring all results 
from ATLAS, CMS and the Tevatron experiments together and test not only the SM but also 
specific models beyond. Such fits by theorists, using various parametrizations, were performed 
previously in 16 43 . Here, we go beyond these works by including all publicly available data 



as of the end of 2012. In particular we take into account the updates presented at the Hadron 
Collider Physics Symposium in Nov. 2012 (HCP2012) [44] and at the Open Session of the 
CERN Council in Dec. 2012 [45]. 

Our parametrization is as follows. We treat the couplings to up-type and down-type 
fermions, Cjj and C*d, as independent parameters (but we only consider the case Cl = Co, and 



1 Although ATLAS finds a lower value of Mh — 123.5 GeV in the ZZ channel, the combined value from 
ATLAS is M H = 125.2 ± 0.3 (stat) ± 0.6 (sys) GeV, in very good agreement with M H = 125.8 ± 0.4 (stat) ± 
0.4 (sys) GeV measured by CMS. We find it reasonable to assume that the lower Mh value from the ATLAS 
H — >• ZZ measurement is due to a statistical fluctuation or unknown systematics. 
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we assume that the Cp are family universal). Moreover, we assume a custodial symmetry in 
employing a single Cw = Cz = Cy in Eq. (JT1) . The structure we are testing thus becomes 



C = g 



cos 6^ ^ I 2m w 2m w 2m w 



H. (2) 



In general, the Ci can take on negative as well as positive values; there is one overall sign ambi- 
guity which we fix by taking Cy > 0. Even in this restricted context, various types of deviations 
of these three Cj from unity are possible in extended theories such as Two-Higgs-Doublet Mod- 
els (2HDMs), models with singlet-doublet mixing, and supersymmetric models such as the 
Minimal Supersymmetric Standard Model (MSSM) and the Next-to-Minimal Supersymmetric 
Standard Model (NMSSM). 

In addition to the tree-level couplings given above, the H has couplings to gg and 77 that 
are first induced at one loop and are completely computable in terms of Cjj, Cp and Cy if 
only loops containing SM particles are present. We define C g and C 7 to be the ratio of these 
couplings so computed to the SM (i.e. C\j = Cd = Cy = 1) values. However, in some of our 
fits we will also allow for additional loop contributions AC S and AC 7 from new particles; in 
this case C g = C g + AC g and C 7 = C 7 + AC 7 . The largest set of independent parameters in 
our fits is thus 

Cu, Cd, Cy, AC 9 , AC 7 . (3) 

In this study, we focus on models in which the Higgs decays only to SM particles, in 
particular not allowing for invisible (e.g. H — y X1X11 where Xx is the lightest SUSY particle) or 
undetected decays (such as H — y aa, where a is a light CP-odd, perhaps singlet scalar). This 
approach, when we allow in the most general case for the Cu, Cd, Cy, C 7 and C g couplings 
to be fully independent, encompasses a very broad range of models, including in particular 
those in which the Higgs sector consists of any number of doublets + singlets, the only proviso 
being the absence of decays of the observed ~ 125 GeV state to non-SM final states. (A fit 
for invisible Higgs decays was performed early on in [46].) This approach however does not 
cover models such as composite models and Higgs-radion mixing models for which the VVH 
coupling has a more complicated tensor structure than that given in Eq. (|2]). Our procedure 
will also be inadequate should the observed signal at ~ 125 GeV actually arise from two or 



more degenerate Higgs bosons (see e.g. 47,481). Although the success of our fits implies that 



there is no need for such extra states, the explicit tests for degenerate states developed in 49 
should be kept in mind as a means to test directly for two or more Higgs bosons contributing 
to the signal at 125-126 GeV. 

This paper is organized as follows. The experimental inputs and our fitting procedure are 
described in Section 2. The results of three generic fits are presented in Section 3 together with 
the results of a fit in Two-Higgs-Doublet models. Section 4 contains our conclusions. 
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2 Experimental inputs and fitting procedure 



We perform fits employing all production/decay channels for which results are available 
from the ATLAS and CMS collaborations at the LHC, as well as the Tevatron CDF+DO Higgs 
results. The experimental results are given in terms of signal strengths /i(A, Y), the ratio of the 
observed rate for some process X — > H — > Y relative to the prediction for the SM Higgs. Often 
it is the case that several production processes contribute to a given experimental channel. 
For example, both vector boson fusion and gluon fusion can contribute to the "VBF" channels 
(or "categories") that are defined by a given set of experimental cuts. In comparing theory 
to experiment it is thus important to incorporate the estimates from the experiments of the 
relative contributions of the theoretically distinct production/decay processes. The values for 
the signal strengths in the various (sub) channels as reported by the experiments and used in 
this analysis, together with the estimated decompositions into production channels are given 
in Tables QM 

We adopt the simple technique of computing the x 2 associated with a given choice of the 
input parameters following the standard definition: 

2 ^ k - ^ 2 (A\ 

where k runs over all the experimentally defined production/decay channels employed, is 
the observed signal strength for channel k, JI k is the value predicted for that channel for a given 
choice of parameters and is the experimental error for that channel. The JI k associated 
with each experimentally defined channel is further decomposed as 



( 5 ) 



where the T % k give the amount of contribution to the experimental channel k coming from the 
theoretically defined channel % and juj is the prediction for that channel for a given choice of Cjj, 
Cd, Cy and (for fits where treated as independent) C 1 and C g . For the computation of the % in- 
cluding NLO corrections we follow the procedure recommended by the LHC Higgs Cross Section 



Working Group in 60 . In particular we include all the available QCD corrections for C„ using 



HIGLU |6TJ[62] and for C 1 using HDECAY 62 63 , and we switch off the electroweak corrections. 
The T k depend on the specific analysis and hence differ from experiment to experiment. Often, 
the T k are determined from simulations of a SM Higgs signal. In some cases, the experiments 
have done the unfolding of theoretical vs. experimental channels from the data and provide di- 
rectly experimental results for the theoretically relevant //j's and the correlations between them. 



With this framework programmed, our fitting procedure is as follows. We first scan over 
a fine grid of the free parameters of the scenario considered, for example, Cu, Co, Cy with 
C g , C 7 = C g , C 1 as computed from the SM-particle loops (this will be Fit II below). We obtain 
the value of \ 2 associated with each point in the grid and thus determine the values of the 
parameters associated with the approximate minimum (or minima). To get the true minimum 
X 2 , Xmin' an d the associated "best-fit" values and the one-standard deviation (la) errors on 
them we employ MINUIT [64]. (The errors on parameters which are not input, i.e. C g and C 7 , 



3 



Channel 


Signal strength p 


M H (GeV) 


Production mode 
ggF VBF VH ttH 


H -> 77 (4.8 fb" 1 at 7 TeV + 13.0 fb" 1 at 8 TeV) < 


1 


/z(ggF + ttH,7 T ) 
/x(VBF + VH, 77 ) 


1.85 ±0.52 
2.01 ± 1.23 


126.6 
126.6 


100% 

60% 40% 


H ^ ZZ (4.6 fb" 1 at 7 TeV + 13.0 fb" 1 at 8 TeV) |6[|ll| 


Inclusive 1.01±g;|g 


125 87% 7% 5% 1% 


H -»■ WW (13.0 fb" 1 at 8 TeV) 8 11 


evfiu 1.42±g;|| 


125.5 95% 3% 2% 


if -»■ 66 (4.7 fb" 1 at 7 TeV + 13.0 fb" 1 at 8 TeV) 11, 


50 


VH tag -0.39 ± 1.02 


125.5 




100% 


H — >■ rr (4.6 fb" 1 at 7 TeV + 13.0 fb -1 at 8 TeV) [51] 


/i(VBF + VH, rr) 


2.41 ± 1.57 
-0.26 ± 1.02 


125 
125 


100% 

60% 40% 



Table 1: ATLAS results as employed in this analysis. The correlations included in the fits are 
p = —0.37 for the 77 and p = —0.50 for the rr channels. 



Channel 


Signal strength p 


M H (GeV) 


Production mode 
ggF VBF VH ttH 


H -> 77 (5.1 fb -1 at 7 TeV + 5.3 fb -1 at 8 TeV) |2[|5||12| 


/i(ggF + ttH, 77 ) 
/i(VBF + VH, 77 ) 


0.95 ±0.65 
3.77 ± 1.75 


125.8 
125.8 


100% 

60% 40% 


H -»■ ZZ (5.1 fb -1 at 7 TeV + 12.2 fb" 1 at 8 TeV) 7 


12| 


Inclusive 0.81^ 


125.8 


87% 7% 5% 1% 


H WW (up to 4.9 fb" 1 at 7 TeV + 12.1 fb" 1 at 8 TeV) 


1( 


), 12,52] 


0/1 jet 
VBF tag 
VH tag 


77+ u - 27 

u - ' ' -0.25 

-0 05 +0 - 74 

u - ud -0.55 
u - oi -1.94 


125.8 
125.8 
125.8 


97% 3% 
17% 83% 

100% 


H -> bb (up to 5.0 fb" 1 at 7 TeV + 12.1 fb" J 


at 8 TeV) [12, 


53 


54 


VH tag 
ttH tag 


1 01 +0.65 
i - oi -0.60 

-O.SO^; 1 " 


125.8 
125.8 




100% 

100% 


H — >• tt (up to 5.0 fb" 1 at 7 TeV + 12.1 fb" 


1 at 8 TeV) (l2, 


55 


56 


0/1 jet 
VBF tag 
VH tag 


85 +abS 
0.82t f 5 


125.8 
125.8 
125.8 


76% 16% 7% 1% 
19% 81% 

100% 



Table 2: CMS results as employed in this analysis. The correlation included for the 77 channel 
is p = -0.54. 
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Channel 


Signal strength \i 


(GeV) 


Production mode 
ggF VBF VH ttH 


H — > 77 


59 




Combined 


6-141^ 


125 


78% 5% 17% 


H^WW 59 




Combined 


U - OJ -0.81 


125 


78% 5% 17% 


H ^bb 


14] 


VH tag 


l.OU_ 73 


125 


100% 



Table 3: Tevatron results for up to 10 fb 1 at ^fs = 1.96 TeV, as employed in this analysis. 



are determined from the grid data.) For plotting distributions of x 2 as a function of any one 
variable, we use the above grid data together with the best fit value, to profile the minimal 
X 2 value with respect to the remaining unconstrained parameters. The 68%, 95% and 99.7% 
confidence level (CL) intervals are then given by x 2 = Xmm + 1' +4 an d +9, respectively. Two- 
dimensional x 2 distributions are obtained analogously from a grid in the two parameters of 
interest, profiling over the other, unseen parameters; in this case, we show contours of x 2 corre- 
sponding to the 68% ( X 2 = xLx + 2-30), 95% ( X 2 = ^ + 6.18) and 99.7% ( X 2 = Xmm + H-83) 
confidence levels for 2 parameters treated jointly. 

Before presenting our results, a couple of comments are in order. First of all, we stress that 
in models of new physics beyond-the-SM (BSM), both the branching ratios and the production 
cross sections and distributions (and indeed the number of Higgs particles) may differ from SM 
expectations. For any BSM interpretation of the Higgs search results it is absolutely crucial to 
have as precise and complete channel-by-channel information as possible [65]. Unfortunately, 
not all the experimental analyses give all the necessary details. Below we comment on how we 
use the currently available information from the experiments. The ideal case would of course 
be that the experiments consistently do the unfolding of theoretical vs. experimental channels 
from the data as mentioned above and always provide directly the experimental results for the 
theoretically relevant /ij's (see Eq. ([5])) and the correlations between them. 



ATLAS 



H — > 77: we take our information from Fig. 4 of [4j. This figure shows the results after 
unfolding to obtain the experimental results for the /x's as defined theoretically. Fig. 4 
does make the approximation that VBF and VH can be lumped together (i.e. have the 
same efficiencies after cuts) and that ggF and ttH can be similarly lumped together (note 
that ttH contributes less than 1%). We fit the 68% CL contour assuming that the A% 2 
follows a bivariate normal distribution]^] With this, the correlation p = —0.37 between 
the ggF and VBF+VH channels is automatically taken into account. We also note that 
while Fig. 4 of (4] is for 126.6 GeV, Fig. 12 (right) in the same paper shows that there is 
a broad "plateau" as a function of the mass when the energy scale uncertainty is taken 
into account, implying that the results should not depend too much on the mass. 



We thank Guillaume Drieu La Rochelle for providing this fit, cf. Table 4 in version 2 of 39 
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j+0.53 
-0.48 



H — > ZZ: the signal strength in this channel reported by ATLAS j6 11 is /x = 1.3 
with a best fit mass of M H = 123.5 ± 0.9 (stat) ± 0.3 (sys) GeV. At M H = 125 GeV, the 
signal strength is /i = l.Ol^o^g, see Fig. 10 in 11 . Assuming that the discrepancy in the 



Higgs mass determined from the 77 and the 4 lepton final states is due to a statistical 
fluctuation (rather than unknown systematics) we use p(ZZ) at Mh = 125 GeV, i.e. 
close to the combined best fit mass from ATLAS, in our fits. 

Alternatively, one could rescale the value of p = 1. 3^048 at Mh = 123.5 GeV for a 
Higgs mass of 125 GeV. This would give p{ZZ) = 1.15±g;^| at M H = 125 GeV (or 
p(ZZ) = l.lllo 4I & t Mff = 125.5 GeV). We checked that taking this alternative approach 
has only marginal influence on our results. 

Regarding the decomposition in production modes, no statement is made in the conference 
note or paper. However, as it is an inclusive analysis, we take the relative ratios of 
production cross sections for a SM Higgs as a reasonable approximation. To this end, we 



use the ratios given by the LHC Higgs Cross Section Working Group 57 . 



H — > WW: we adopt relative contributions of 95% ggF and 5% VBF |8|. We do not 
include any result for 7 TeV because the update presented at HCP is a combination of 7 
and 8 TeV. 

H — » rr: ATLAS provides only an overall signal strength with no information on the 
decomposition with respect to production modes. However, the conference note 51 con- 
tains the results of unfolding to the theory-level p's via a plot (Fig. 19) of the experimental 
results for p gg p x B/Bsu versus Pvbf+vh x B/B^m at Mjj = 125 GeV. We fit the con- 
tours with the same procedure as for H — > 77. The correlation (included in the fit) is 
p = -0.50. 



CMS 



H — > 77: we follow the same procedure as for ATLAS H — >■ 77, using Figure 11 from 12 . 
The correlation (included in the fit) is p = —0.54. 

H — > ZZ: no decomposition with respect to production modes is given in the conference 
note or paper. As it is a fully inclusive analysis, we use the relative ratios of production 
cross sections given by the LHC Higgs Cross Section Working Group [57] as a good 



approximation 58 . 



H — > WW: the information provided in the conference note and papers is incomplete; 
our decomposition into production modes is based on [58). Our combination (weighted 
mean) agrees within 9% with that given by CMS (/i CO mb = 0.64±0.24 instead of 0.70^ 



K>.24^ 
-0.23 J 



H — > bb: as there is no information on possible contaminations by other production 
modes, we assume 100% VH or 100% ttH production for the respective categories. 

H — > tt: for the 0/1 jet and VBF tag categories we extract the decomposition into 
production modes from |56|, assuming that there is no significant change in the efficiencies 
between Mh = 125 GeV and Mh = 125.8 GeV. We use the efficiencies from the first 
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three categories (fir h + X, er h + X and e/i + X) because they are the most sensitive ones; 
they lead to very similar decompositions which we use in our analysis. Our combination 
(weighted mean) agrees within 6% with that given by CMS (/i CO mb = 0.83 ± 0.49 instead 
of 0.88j£g). 



Tevatron 



H — > 77 and H — > WW: no decomposition into production modes is given by the 
experiments. We assume that the analyses are inclusive and we thus employ the ratios of 
the theoretical predictions for the (SM) Higgs production cross sections. 



H — > bb: we use the new results from HCP2012 14 assuming 100% VH. 



3 Results 



3.1 General coupling fits 

Fit I: C v = C D = C v = 1, AC g and AC 7 free 

For a first test of the SM nature of the observed Higgs boson, we take Cy = Cy, = Cy — 1 
{i.e. quark, lepton and W, Z vector boson couplings to the Higgs are required to be SM-like) but 
we allow for additional new physics contributions to the 77 and gg couplings, parameterized by 
AC g and AC 7 , coming from loops involving non-SM particles or from anomalies. This fit, which 
we refer to as Fit I, is designed to determine if the case where all tree-level Higgs couplings are 
equal to their SM values can be consistent with the data. For example, such a fit is relevant in 



the context of UED models where the tree-level couplings of the Higgs are SM-like 66,67 . 

Figure [T] displays the results of this fit in the AC g versus AC 7 plane. The best fit is obtained 
for AC 7 ~ 0.43, AC g ~ —0.09, and has xLm = 12.31 for 19 degrees of freedom (d.o.f.), giving 
a p- value of 0.87. The results of this fit are summarized in Table |4| together with the results 
of the other fits of this section. 

We note that the SM (i.e. Cy = C D = Cy = 1, AC g = AC 7 = 0) has x 2 = 20.2 and is 
hence more than 2a away from the best fit in Fig. [Tj The number of degrees of freedom for the 
SM fit is 21, implying a p-value of 0.51. The largest x 2 contributions come from the H — > 77 
channels from ATLAS (A X 2 = 5.06), CMS (A X 2 = 3.36) and Tevatron (A^ 2 = 2.60), followed 
by the VBF result for H -> WW from CMS with A X 2 = 2.01. 



Fit II: varying C v , C D and C v (AC 7 = AC g = 0) 

Next, we let CV, Cp, Cy vary, assuming there are no new particles contributing to the 
effective Higgs couplings to gluons and photons, i.e. we take AC 7 = AC g = implying C g = C g , 
C 7 = C 7 as computed from the SM-particle loops. The results for the one-dimensional and 
two-dimensional x 2 distributions are shown in Figs. [2] and [3] The value of Cy is rather well 
determined to be close to unity. It is intriguing that the best fit of Cy is indeed just slightly 
below 1, as any model with only Higgs doublets or singlets requires Cy < 1. The best fit values 
for Co and C\j are SM-like in that they have magnitudes that are close to one. However, the 
best fit Cy value is opposite in sign to the SM Higgs case. The preference for Cy < is at 
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1 



0.5- 




-1 1 1 1 ' 1 

-1 -0.5 0.5 1 

AC 7 

Figure 1: Two parameter fit of AC 7 and AC g , assuming Cjj = Cd = Cy = 1 (Fit I). The red, 
orange and yellow ellipses show the 68%, 95% and 99.7% CL regions, respectively. The white 
star marks the best-fit point AC 7 = 0.426, AC g = -0.086. 

the level of 2.6a — see the first plot in Fig. [2j This results from the fact that an enhanced 77 
rate (as observed in the experimental data) is obtained by changing the sign of the top-loop 
contribution so that it adds, rather than subtracts, from the W^-loop. In contrast, in the case of 
Cd almost equally good minima are found with Cd < and Cd > 0. Details on the minima in 
different sectors of the (Cu, Cd) plane are given in Table[5] Note that, for the best fit point, the 
resulting C 7 and C g are in good agreement with the result of Fit I above, for which C 7 = 1.43 
and C g = 0.91. Here, however, the enhanced C 1 value derives from Cu < rather than from 
AC 7 7^ 0. The best fit results are again tabulated in Table |4j 

A negative sign of C\j — while maintaining a positive sign of m t — is actually not easy to 
achieve. (A sign change of both C\j and m t would have no impact on the top quark induced 
loop amplitudes.) It would require that m t is induced dominantly by the vev of a Higgs boson 
which is not the Higgs boson considered here. Hence, we have Cu > in most models, implying 
that it is important to study the impact of this constraint on our fits. The fit results when 
requiring Cu, Cd > are shown in the left two plots of Fig. [4] and the top row of Fig. |5j see 
also Table [5j We observe that for this quadrant the results are consistent with SM expectations 
(i.e. within ~ la). Interestingly the fit is not better than the SM itself: Xmm = 18-66 for 
21 — 3 = 18 d.o.f., corresponding to p — 0.41. 

Another possible model constraint is to require Cy < 1 (recall that Cy > by convention). 
This constraint applies to any model containing only Higgs doublets and singlets. The Id 
results for the combined requirement of Cu, Cd > and Cy < 1 are shown in the right two 
plots of Fig. |4j and in the bottom-row plots of Fig. [5| We observe that the best fit values for 
Cu and Cd are only slightly shifted relative those found without constraining Cy, and that 
accordingly the C 7 = C 7 and C g = C g at the best fit point are only slightly shifted. However, 
the Cy < 1 constraint does severely change the upper bound on C 7 , which for Cu > and 
AC 7 = mostly depends on the W-boson loop contribution. The apparent sharpness of the 
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Figure 3: Two-dimensional y 2 distributions for the three parameter fit, Fit II, of Cjj, Cr>, Cy 
with C 7 = C 7 and C g = C g as computed in terms of Cjj, Co, Cy. The red, orange and yellow 
ellipses show the 68%, 95% and 99.7% CL regions, respectively. The white star marks the 
best-fit point. Details on the minima in different sectors of the (CV, Co) plane can be found in 
Table B 
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Cu,d,v C lt g Cjj,d,V C 7 , 9 



Figure 4: One-dimensional y 2 distributions for the three parameter fit, Fit II, but imposing 
Cjj > 0, Cd > 0; the left two plots allow for Cy > 1 (Xmin = 18.66), while in the right two 
plots C v < 1 (xLn = 18-89). 




Figure 5: Two-dimensional x 2 distributions for the three parameter fit, Fit II, as in Fig. [3] but 
with Cjj > 0, Cd > 0, Cy > 0. The upper row of plots allows for Cy > 1, while in the lower 
row of plots Cy < 1 is imposed. 



10 



Fit 


I 


II 


III, 1st min. 


III, 2nd min. 


Cu 


1 


-0.86^ 


-0.06 ± 1.30 


0.06 ±1.30 


c D 


1 




1 oo +a28 

l.UU_ 2 Q 


-1 oo +0 - 26 

l.UU__g 28 


C v 


1 


qs +0 ' 12 


Q3+ 012 

U - JO -0.14 


Q3+ 012 

U - JO -0.14 


AC 1 


n 43+ - 17 


— 


1 6 +0 - 38 

U ' iu -0.36 


91 + - 37 
u - zi -0.39 




-0.09 ±0.10 


— 


0.83^ 


0.83i?# 




1 4Q+ - 17 


1.43 ±0.17 


1 QC+0.26 
±.OU_g 23 


1 QC+0.26 
!.OU_g 23 


c 9 


0.91 ±0.10 


Q2 +a17 


u.yj_ 23 


Q c i +a26 
u.yj_ 23 


A mm 


12.31 


11.95 


11.46 


11.46 


X min /d.o.f. 


0.65 


0.66 


0.72 


0.72 



Table 4: Summary of results for Fits I— III. For Fit II, the tabulated results are from the best 
fit, cf. column 1 of Table |5j 



Sector 


C V < 0, C D > 


Cu, C D < 


Cu, C D > 


Cu 


-0.86182$ 


-0.91^ 


u -° d -0.13 


c D 


0.99^;1 


-0.98^; 2 ? 


o s^ - 22 

U.OO_ 2i 


C v 


Q c i+ a12 


Q4+ - 12 
u - J ^-0.13 


i nfi+o.n 

i - uu -0.12 


Cry 


1.43 ±0.17 


1 AQ+0.16 
L - w -0.17 


i.n±S:ii 


c 9 


Q2 +ai7 


qi+ a17 

u - yi -0.15 


o.85±g:i| 


v 2 . 
Amin 


11.95 


12.06 


18.66 


Xmin/ d - - f - 


0.66 


0.67 


1.04 



Table 5: Results for Fit II in different sectors of the (Cu, Co) plane. 



boundary in the C g vs. C 7 plane is a result of the fact that these two quantities really only 
depend on Cu for Cy = 1. 

Finally note that it has been shown in 



68 , 69 that single top production in association 



with a Higgs is greatly enhanced when Cu,Cy have opposite signs. Thus, the possibility of 
Cu < should be further scrutinized by precision measurements of the single top production 
cross section at the LHC. 



Fit III: varying C v , C D , C v , AC 7 and AC g 

Finally, in Fit III, we allow the AC 9 and AC 7 additions to C g and C 7 , fitting therefore to 
five free parameters: Cu, Co, Cy, AC g , and AC 7 . The associated Id and 2d plots are given 
in Figs. [6] and [7} There are two main differences as compared to Fit II. On the one hand, the 
preference for C 7 > 1 does does not necessarily imply a negative value for Cu, since a positive 
value for AC 7 can contribute to tin iiicrGciSG in Cry even when the top-quark loop interferes 
destructively with the W loop. (This is obviously already expected from Fit I.) On the other 
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hand, both Cu and AC g feed into the effective C g , and if one of them is large the other one 
has to be small to result in a near SM-like gg — > H cross section. This anti-correlation between 
\Cu\ and AC g can be seen in the center-top plot in Fig. [7j The best fit is actually obtained for 
Cu ~ 0, with AC g ~ 1 in order to compensate for the very suppressed top- loop contribution 
to ggF. However, it is also apparent that the minimum at Cu = is quite shallow (cf. the top 
left plot in Fig. [6]) and that a fit with Cu ~ 1 with small AC g is well within the 68% contour 
(as should indeed be the case for consistency with Fits I and II). 

We also note that at the best fit, i.e. that with Cu ~ 0, one finds C 7 ~ C 7 > 1 by 
virtue of the fact that the W loop is not partially cancelled by the top loop and only a small 
AC 7 ~ 0.16-0.21 is needed to further enhance the 77 final state and bring 7/(77 ) into agreement 
with observations; see top-right and bottom-right plots of Fig. [7j If we move to the SM value 
of Cu = 1 then AC 7 ~ 0.45 is needed to fit the 77 rate. The best fit results are tabulated in 
Table S 

A way to lift the degeneracy in Cu and AC g would be to have an independent determination 
of Cu- This might be achieved by an accurate measurement of the ttH channel, as illustrated 
in Fig. [8j This figure assumes that //(ttH) will eventually be measured with 30% accuracy - 
more concretely, the figure assumes //(ttH) = 1 ± 0.3. This is certainly a very challenging task. 



For comparison, CMS currently gives //(ttH) « — 0.8_ 1 ; 8 12 . Finally, as mentioned above, Cu 



may also be constrained by the associated production of a single top and a Higgs 68 , 69 
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El 

-2-1 1 2-2-1 1 2 

Figure 8: Lifting of the degeneracy in Cjj and AC 9 in Fit III when tiH is measured to 30% 
accuracy (/x(ttH) = 1 ±0.3). These two plots should be compared to the top left and top middle 
plots of Fig. [7j See text for details. 
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Type I and II 


Type I 


Type II 


Higgs 


VV 


up quarks 


down quarks & 
leptons 


up quarks 


down quarks & 
leptons 


h 


sin(/3 — a) 


cos a / sin (3 


cos a/ sin/3 


cos a/ sin /3 


—sin a/ cos/3 


H 


cos(/3 — a) 


sin a j sin 


sin a / sin /? 


sin a/ sin /3 


cos a / cos /3 


A 





cot (3 


— cot 


cot (3 


tan /3 



Table 6: Tree-level vector boson couplings Cy (V = W, Z) and fermionic couplings C F i nor- 
malized to their SM values for the Type I and Type II two-Higgs-doublet models. 

3.2 Two-Higgs-Doublet Model 

So far our fits have been model- independent, relying only on the Lagrangian structure of 
Eq. g. Let us now turn to the concrete examples of Two-Higgs-Doublet Models (2HDMs) of 
Type I and Type II. In both cases, the basic parameters describing the coupling of either the 
light h or heavy H CP-even Higgs boson are only two: a (the CP-even Higgs mixing angle) 
and tan/3 = v u /vd, where v u and Vd are the vacuum expectation values of the Higgs field that 
couples to up-type quarks and down-type quarks, respectively. The Type I and Type II models 
are distinguished by the pattern of their fermionic couplings as given in Table [6j The SM limit 
for the h (H) in the case of both Type I and Type II models corresponds to a = (3 — %/2 
(a = (3). In our discussion below, we implicitly assume that there are no contributions from 
non-SM particles to the loop diagrams for C 7 and C g . In particular, this means our results 
correspond to the case where the charged Higgs boson, whose loop might contribute to C 7 , is 
heavy. 

The results of the 2HDM fits are shown in Fig. [9] for the case that the state near 125 GeV 
is the lighter CP-even h. The figure also applies for the case of the heavier H being identified 
with the ~ 125 GeV state with the replacement rules given in the figure caption]^] Note that 
the convention Cy > implies sin(/3 — a) > for the h and cos(/3 — a) > for the H. Moreover, 
the requirement tan/3 > restricts (3 G [0, tt/2]. The best fit values and la ranges for a and (3, 
together with the corresponding values for CV, Cd, Cy, C g and C g , are listed in Table [7| These 
numbers are again for the case of h being the state near 125 GeV. Replacing h by H amounts 
to a shift in a -> a + vr/2; thus we find a = 6.07±g;g| (cos a = 0.98 ± 0.02) for the 2HDM-I and 
a = 6.14+°;^ (cos a = 0.99i°;^) for the 2HDM-II, while the values for tan/3, C v , C D , C v , etc. 
do not change. 

Note that for both the Type I and the Type II model, the best fits are quite far from the SM 
limit in parameter space. In particular, since we do not include any extra loop contributions 
to C 7 , we end up with negative Cy close to —1 as in Fit II. Demanding Cy > (i.e. cos a > 
for h, sin a > for H), one ends up in a long 'valley' along the decoupling limit where the 
Higgs couplings are SM like, see Fig. [9j this is however always more than 2a away from the 
best fit. Furthermore, solutions with very small tan/3 < 1 are preferred at more than 2a. Since 
such small values of tan/3 are rather problematic (in particular tan/3 < 0.5 is problematical 

3 Since the ^125 GeV state clearly couples to WW, ZZ we do not consider the case where the A is the only 
state at ~ 125 GeV. We also do not consider the cases where the ~ 125 GeV peak comprises degenerate (h, H), 
(h, A) or (H, A) pairs. 
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cos a cos a 



Figure 9: 2HDM fits for the h in the Type I (left) and Type II (right) models. The upper row 
shows the fit results in the (3 — tt/2 vs. a plane, while the lower row shows the sin/3 vs. cos a 
plane. The dashed lines indicate the SM limit. The same results are obtained for the heavier 
H with the replacements (3 — tt/2 — > (3 and a — > a + tt/2 (sin (3 — > — cos (3, cos a — > sin a). 



for maintaining a perturbative magnitude for the top-quark Yukawa coupling) we also give in 
Table [7] the corresponding fit results requiring tan > 1. These results come quite close to the 
SM limit, and accordingly have a xLin °f about 19-20 (recall that for the SM we find x 2 — 20.2). 
2HDMs with tan (3 > 1 hence do not provide a better fit than the SM itself. 

A couple of more comments are in order. First, an important question that we leave for 
future work is whether other - - e.g. stability, unitarity, perturbativity (SUP) and precision 
electroweak (PEW) — constraints are obeyed at the best-fit points, or the 68% CL regions. 



Here we just note that according to Fig. 1 of 70 , the SUP and PEW constraints do not seem 



problematic for Type II, but may play a role for Type I models at low tan (3. 

Second, the best fits correspond to very small tan/3 (small (3) values that are potentially 
constrained by limits from B-physics, in particular from AM# S and Z — >■ bb . The B-physics 
constraints are summarized in Figs. 15 and 18 of [7l] for Type II and Type I, respectively. 
Figure 18 for Type I places a lower bound on tan (3 as a function of the charged Higgs mass 
which excludes small tan (3 < 1 unless the charged Higgs is very heavy, something that is 
possible but somewhat unnatural. Figure 15 for Type II places a substantial lower bound on 
the charged Higgs mass for all tan (3, but such a constraint does not exclude the 68% CL region. 
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Fit 


2HDM-I 


2HDM-II 


2HDM-I, tan p > 1 


2HDM-II, tan > 1 


a [rad] 


4 c,n +0 - 09 

^• ou -0.08 


4 c;fi+o.i5 

^• oo -0.14 


r 07+l.H 


6-28l8:s3 


P [rad] 


u - z ^-0.10 


17+ - 12 
V - L ' -0.17 


[tt/4, tt/2] 


1 ^Pl +0m 

1.<JU_Q 7 g 


cos a 


n 91 +0.09 
u - zi -0.08 


—0 1 c;+ - 15 

u - io -0.13 


0.61185? 


l-00_ .67 


tan P 


24 +0 ' 08 

u - z ^-0.10 




[1, +oo[ 


[1, +oo[ 


Cu 


U - JU -0.19 


-o.87j5:lS 


87+ ' 17 

u -°' -0.15 


1 02 +0 ' 05 

± - uz -0.07 


c D 


-0 qo +a17 

U - JU -0.19 


1.00_ .oi 


87 +0 - 17 

u -°'-0.15 


n 04+0.13 


C v 


0.90 ±0.07 


u - yd -0.12 


n qq+ - 01 

u - JJ -0.04 


1.00_ .05 


Cry 


-1 07+0.09 
- 1 -- ' -0.10 


1.4418;?* 


1.03_o.o6 


1 m+ - 01 

i - ui -0.09 


c 9 


0.90l°;^ 


0.921°;}? 


0.8718:11 


n qq+o.08 

u - yy -0.04 


Amin 


12.20 


11.95 


19.43 


19.88 



Table 7: Summary of fit results for the h in 2HDMs of Type I and Type II. 



Third, we remind the reader that in the 2HDMs, the soft ^-symmetry-breaking m 2 2 and 



the other Higgs masses (M h , M H and M4) are independent parameters. It is thus possible to 
have either or Mh ~ 125 GeV without violating constraints from direct searches for the 
charged Higgs whose mass is related to uia- However, in the case of Mh ~ 125 GeV, one has 
to avoid the LEP limits for the lighter h, which severely constrain the h coupling to ZZ in case 

So either M h > 114 GeV for M H « 125 GeV, or sin 2 (/3 - a) needs to 
100 GeV, or sin 2 (/3 - a) < 0.1 for M h < 90 GeV). 
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of M h < 114 GeV 

be small (e.g. sin 2 (/3 — a) < 0.3 for M h 
The Ax 2 distributions of sin 2 (/3 — a) for Type I and Type II with Mh ~ 125 GeV are shown 
in Fig. [lOj Interestingly, around the best fit the h coupling to ZZ is sufficiently suppressed to 
allow for Mh of the order of 100 GeV (or lower in Type II). 





sin 2 (/? - a) 



0.2 0.4 0.6 0.8 

sin {(3 — a) 



Figure 10: A% 2 distribution of sin 2 (/3 — a) in the Type I (left) and Type II (right) models for 
the case that H is the observed state near 125 GeV. 



16 



4 Summary and Conclusions 



We assessed to which extend the currently available data on the Higgs-like scalar constrain 
the Higgs couplings. To this end we performed fits to all public data from the LHC and the 
Tevatron experiments. 

First, we employed a general parametrization of the Higgs couplings based on a SM-like 
Lagrangian, but allowing for extra contributions to the loop-induced couplings of the Higgs-like 
scalar to gluons and photons. While the SM does not provide a bad fit (x 2 /d.o.f. = 0.96), it is 
more than 2a away from our best fit solutions. The main pull comes from the enhanced H — > 77 
rates observed by ATLAS and CMS, as well as from the Tevatron experiments. The best fits 
are thus obtained when either Cu ~ —1 {i.e. opposite in sign to the SM expectation) or there is 
a large BSM contribution to the 77 coupling of the Higgs. In short, significant deviations from 
the SM values are preferred by the currently available data and should certainly be considered 
viable. Since having Cu ~ —1 (in the Cy > convention) is not easy to achieve in a realistic 
model context, and leads to unitarity violation in WW — > ti scattering at scales that can be as 



low as 5 TeV 73 , 74 , it would seem that new physics contributions to the effective couplings 
of the Higgs to gluons and photons are the preferred option. (The possibility of a second, 
degenerate Higgs boson contributing to the observed signal remains another interesting option, 
not considered here.) 

Second, we examined how well 2HDM models of Type I and Type II fit the data. We found 
that it is possible to obtain a good fit in these models with sin(/9 — a) (cos(/3 — a)), in the 
h (H) cases, respectively, not far from 1. However, the best fit values for the individual Cu, 
Co-, C 7 and C g parameters lie far from their SM values. Further, the best fits give tan j3 < 1, 
which is disfavored from the theoretical point of view if we want perturbativity up to the GUT 
scale. Requiring tan/3 > 1 (or simply Cu > 0) pushes the fit into the SM 'valley' and no 
improvement over the pure SM solution is obtained. In particular the x 2 obtained in this 
region is substantially larger than that for the best fit, and not far from the \ 2 found for the 
SM. 

We once again refer the reader to Tables [4j [5] and [7] which summarize the best fit values 
and la errors for the parameters for the various cases considered. In Fig. 11 we show some 
of these results graphically. Moreover, in order to assess the physics associated with our best 
fit points, we give in Tables [8] and [9] the values of the derived (theory level) signal strengths 
ju(ggF,77), £(ggF,ZZ), /2(ggF,66), /2(VBF,7 7 ), j2(VBF,ZZ), and /2(VBF, 66) for the best fit 
point in the various fits we have considered. (These are a complete set since for the models 
we consider $(X,tt) = /2(X,66), fi(X, WW) = j2(X,ZZ) and yu(VBF, Y) = /2(VH,F).) We 
see that in the general case both /2(ggF,77) and /2(VBF, 77) are enhanced by factors 1.7-2.1 
(1.8-1.9 in 2HDMs), while the other signal strengths tend to be < 1. When demanding Cu > 
without allowing for extra contributions from new particles, then only very small enhancements 
of /i(VBF,77) and /2(VBF, ZZ) of the order of 1.2-1.3 are found. 

Last but not least, we strongly encourage the experimental collaborations to make as com- 
plete as possible channel-by-channel information (including the important decomposition into 
production modes) available, in order to allow for reliable tests of non-standard Higgs scenarios. 
The information currently given by ATLAS and CMS for the 77 signal is an example of good 
practice and should become the standard for the presentation of results for all channels. This 
would immensely help interpretation efforts such as attempted in this paper. The ideal case 
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Figure 11: Graphical representation of the best fit values for Cu, Co, Cy, AC 7 and AC g of 
Table HI The labels refer to the fits discussed in the text. The dashed lines indicate the SM 
value for the given quantity. The x 's indicate cases where the parameter in question was fixed 
to its SM value. 



Fit 


I 


II, C v < 


II, Cu- > 


III 


£(ggF, 77) 

//(ggF, 66) 
£(VBF, 77 ) 
//(VBF, ZZ) 
ju(VBF, 66) 


-i 71+0.33 

L - ' x -0.32 

0.84l°f 7 

0.841™ 
2 p i +0 - 54 

Z.UJ_g 44 

1.00 ±0.02 
1.00 ±0.02 


1 01+0.43 

l-8l_0.41 

0.79 ±0.15 
0.87l°j 7 

1 qo+0.78 

0.92 ±0.30 


1.07 ±0.18 
0.97 ±0.20 

63 +0 - 36 

1 fifi +0.70 
i - oo -0.63 

1 c-Q+0.50 

i.ou_ 046 
0.98 ±0.32 


1 7O+0.36 
L - ' J -0.34 

0.841™ 
0-96±g;g 

1.7^ 

0.82j£fg 
0.931°;! 



Table 8: Summary of ju results for Fits I III. For Fit II, the tabulated results are for the best 
fit with Cu < 0, column 1 of Table [5j and for the case Cu, Cd > 0, column 3 of Table [5j 



Fit 


2HDM-I 


2HDM-II 


2HDM-I, tan > 1 


2HDM-II, tan > 1 


£(ggF, 77) 

/2(ggF, 66) 
/i(VBF, 77 ) 
£(VBF, ZZ) 
/2(VBF, 66) 


1 O/5+0.41 
1.00-0.38 

0.81 ±0.11 

0.80l°;f 7 
1 07+0.77 

J- 8 ' -0.65 
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1 01 +0.41 
i - yi -0.40 

0.79l°;g 
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1.91183" 
0.84l°;i 
0.93l°i° 


n Q7+ - 03 

u - J '-0.09 

Q1+ - 10 
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Table 9: Summary of /1 results for the interpretation in 2HDM models. 



would of course be if the full likelihood distributions were made available. 

18 



Acknowledgements 



We are indebted to Guillaume Drieu La Rochelle for help in extracting the signal strength 
information in the ATLAS 77 and rr and the CMS 77 channels, and to Albert de Roeck 
for detailed discussions on the CMS results. Moreover, we gratefully acknowledge discussions 
with Ada Farilla, Jonathan Hays, Marumi Kado, Tom LeCompte, Koji Nakamura and Mayda 
Velasco. 

This work was supported in part by US DOE grant DE-FG03-91ER40674 and by IN2P3 
under contract PICS FR-USA No. 5872. UE acknowledges partial support from the French 
ANR LFV-CPV-LHC, ANR STR-COSMO and the European Union FP7 ITN INVISIBLES 
(Marie Curie Actions, PITN-GA-2011-289442). GB, UE, JFG and SK acknowledge the hospi- 
tality and the inspiring working atmosphere of the Aspen Center for Physics which is supported 
by the National Science Foundation Grant No. PHY-1066293. GB, JFG and SK also thank 
the Galileo Galilei Institute for Theoretical Physics for hospitality and the INFN for partial 
support. JFG acknowledges the generous hospitality while completing this work of the Kavli 
Institute for Theoretical Physics which is supported by the National Science Foundation under 
Grant No. NSF PHY11-25915. 

References 

[1] G. Aad et al. [ATLAS Collaboration], "Observation of a new particle in the search for the 
Standard Model Higgs boson with the ATLAS detector at the LHC," Phys. Lett. B 716 
(2012) 1 I arXiv: 1207.7214] [hep-ex]]. 

[2] S. Chatrchyan et al. [CMS Collaboration], "Observation of a new boson at a mass of 125 
GeV with the CMS experiment at the LHC," Phys. Lett. B 716 (2012) 30 |arXiv:1207.723~5 
[hep-ex]]. 

[3] Tevatron New Physics Higgs Working Group and CDF and DO Collaborations, "Updated 
Combination of CDF and DO Searches for Standard Model Higgs Boson Production with 
up to 10.0 fb- 1 of Data," |arXiv: 1207.0449] [hep-ex]. 

[4] ATLAS-CONF-2012-168, "Observation and study of the Higgs boson candidate in the two 
photon decay channel with the ATLAS detector at the LHC" (Dec 2012). 

[5] CMS-PAS-HIG-12-015, "Evidence for a new state decaying into two photons in the search 
for the standard model Higgs boson in pp collisions" (July 2012). 

[6] ATLAS-CONF-2012-169, "Updated results and measurements of properties of the new 
Higgs-like particle in the four lepton decay channel with the ATLAS detector" (Dec 2012). 

[7] CMS-PAS-HIG-12-041, "Updated results on the new boson discovered in the search for 
the standard model Higgs boson in the ZZ to 4 leptons channel in pp collisions at y/s = 7 
and 8 TeV" (Nov 2012). 

[8] ATLAS-CONF-201 2-158, "Update of the H WW& evfw Analysis with 13 fb 1 or 
^ = 8 TeV Data Collected with the ATLAS Detector" (Nov 2012). 



19 



[9] ATLAS-CONF-2012-162, "Updated ATLAS results on the signal strength of the Higgs-like 
boson for decays into WW and heavy fermion final states" (Nov 2012). 

[10] CMS-PAS-HIG-12-042, "Evidence for a particle decaying to W + W~ in the fully leptonic 
final state in a standard model Higgs boson search in pp collisions at the LHC" (Nov 2012). 

[11] ATLAS-CONF-2012-170, "An update of combined measurements of the new Higgs-like 
boson with high mass resolution channels" (Dec 2012). 

[12] CMS-PAS-HIG-12-045, "Combination of standard model Higgs boson searches and mea- 
surements of the properties of the new boson with a mass near 125 GeV" (Nov 2012). 

[13] T. Aaltonen et al. [CDF and DO Collaborations], "Evidence for a particle produced in 
association with weak bosons and decaying to a bottom-antibottom quark pair in Higgs 
boson searches at the Tevatron," Phys. Rev. Lett. 109 (2012) 071804 | jarXiv: 1207.6436] 
[hep-ex]]. 

[14] Yuji Enari, U H bb from Tevatron", talk at HCP2012, 14 Nov 2012, Kyoto, Japan, 
http : //kds . kek . jp/ conf erenceDisplay . py?conf Id=10808. 

[15] ATLAS-CONF-2012-127, "Coupling properties of the new Higgs-like boson observed with 
the ATLAS detector at the LHC" (Sep 2012). 

[16] D. Carmi, A. Falkowski, E. Kuflik and T. Volansky, "Interpreting LHC Higgs Results from 
Natural New Physics Perspective," JHEP 1207 (2012) 136 |arXiv:1202.3~l~44] [hep-ph]]. 

[17] A. Azatov, R. Contino and J. Galloway, "Model-Independent Bounds on a Light Higgs," 
JHEP 1204 (2012) 127 |arXiv:1202.3~415] [hep-ph]]. 

[18] J. R. Espinosa, C. Grojean, M. Muhlleitner and M. Trott, "Fingerprinting Higgs Suspects 
at the LHC," JHEP 1205 (2012) 097 larXiv:1202.3~69~7] [hep-ph]]. 

[19] M. Klute, R. Lafaye, T. Plehn, M. Rauch and D. Zerwas, "Measuring Higgs Couplings 
from LHC Data," Phys. Rev. Lett. 109 (2012) 101801 |arXiv:1205.2"699| [hep-ph]]. 

[20] A. Azatov, S. Chang, N. Craig and J. Galloway, "Higgs fits preference for suppressed 
down-type couplings: Implications for supersymmetry," Phys. Rev. D 86 (2012) 075033 
jarXiv: 1206.1058] [hep-ph]]. 

[21] D. Carmi, A. Falkowski, E. Kuflik and T. Volansky, "Interpreting the Higgs," 
|arXiv: 1206.4201] [hep-ph] . 

[22] I. Low, J. Lykken and G. Shaughnessy, "Have We Observed the Higgs (Imposter)?," 
|arXiv: 1207. 1093] [hep-ph] . 

[23] T. Corbett, O. J. P. Eboli, J. Gonzalez-Fraile and M. C. Gonzalez-Garcia, "Constraining 
anomalous Higgs interactions," Phys. Rev. D 86 (2012) 075013 |arXiv:1207.1344| [hep-ph]]. 

[24] P. P. Giardino, K. Kannike, M. Raidal and A. Strumia, "Is the resonance at 125 GeV the 
Higgs boson?," |arXiv: 1207. 1347] [hep-ph]. 



20 



[25] J. Ellis and T. You, "Glob al Analys is of the Higgs Candidate with Mass 125 GeV," JHEP 
1209 (2012) 123 |arXiv: 1207. 1693] [hep-ph]]. 

[26] M. Montull and F. Riva, "Higgs discovery: the beginning or the end of natural EWSB?," 
JHEP 1211 (2012) 018 |arXiv: 1207. 1716 [hep-ph]]. 



[27] J. R. Espinosa, C. Grojean, M. Muhlleitner and M. Trott, "First Glimpses at Higgs' face," 
|arXiv: 1207.1717] [hep-ph] . 

[28] D. Carmi, A. Falkowski, E. Kuflik, T. Volansky and J. Zupan, "Higgs After the Discovery: 
A Status Report," JHEP 1210 (2012) 196 [arXiv:1207.171~8| [hep-ph]]. 

[29] S. Banerjee, S. Mukhopadhyay and B. Mukhopadhyaya, "New Higgs interactions and re- 
cent data from the LHC and the Tevatron," JHEP 1210 (2012) 062 |arXiv:1207.3~58~8| 
[hep-ph]]. 

[30] F. Bonnet, T. Ota, M. Rauch and W. Winter, "Interpretation of precision tests in the 
Higgs sector in terms of physics beyond the Standard Model," |arXiv: 1207.45 99 [hep-ph]. 

[31] T. Plehn and M. Rauch, "Higgs Couplings after the Discovery," Europhys. Lett. 100 (2012) 
11002 (arXiv": 1207.6108| [hep-ph]]. 



[32] J. R. Espinosa, C. Grojean, V. Sanz and M. Trott, "NSUSY fits," |arXiv:1207.7355| [hep- 
ph]. 

[33] D. Elander and M. Piai, "The decay constant of the holographic techni-dilaton and the 



125 GeV boson," arXiv:1208.0546 [hep-ph]. 



[34] A. Djouadi, "Precision Higgs coupling measurements at the LHC through ratios of pro- 
duction cross sections," |arXiv: 1208.3 436 [hep-ph]. 

[35] W. Altmannshofer, S. Gori and G. D. Kribs, "A Minimal Flavor Violating 2HDM at the 
LHC," arXiv: 1210.2465| [hep-ph]. 

[36] B. A. Dobrescu and J. D. Lykken, "Coupling spans of the Higgs-like boson," 
|arXiv: 1210.33421 [hep-ph] . 

[37] S. Chang, S. K. Kang, J. -P. Lee, K. Y. Lee, S. C. Park and J. Song, "Comprehensive 
study of two Higgs doublet model in light of the new boson with mass around 125 GeV," 
|arXiv: 1210.3439] [hep-ph] . 



[38] G. Moreau, "Constraining extra-fermion(s) from the Higgs boson data," arXiv:1210.3977 
[hep-ph] . 

[39] G. Cacciapaglia, A. Deandrea, G. D. La Rochelle and J. -B. Flament, "Higgs couplings 
beyond the Standard Model," [arXiv:1210.8120| [hep-ph]. 

[40] P. Bechtle, S. Heinemeyer, O. Stal, T. Stefaniak, G. Weiglein and L. Zeune, "MSSM 



Interpretations of the LHC Discovery: Light or Heavy Higgs?," arXiv:1211.1955 [hep-ph] 



21 



T. Corbett, O. J. P. Eboli, J. Gonzalez-Fraile and M. C. Gonzalez-Garcia, "Robust Deter- 



mination of the Higgs Couplings: Power to the Data," |arXiv:1211.4580 [hep-ph]. 



E. Masso and V. Sanz, "Limits on Anomalous Couplings of the Higgs to Electroweak 
Gauge Bosons from LEP and LHC," |arXiv:1211.1320"| [hep-ph]. 

A. Azatov and J. Galloway, "Electroweak Symmetry Breaking and the Higgs Boson: Con- 
fronting Theories at Colliders," arXiv: 1212. 1380 [hep-ph]. 



Hadron Collider Physics Symposium 2012 (HCP 2012), 12-16 Nov 2012, Kyoto, Japan, 
http : //www. icepp 7s .u-tokyo . ac . jp/hcp2012/| 

Status of the LHC and Experiments, CERN seminars, 13 Dec 2012, 
|http : / / indico . cern . ch/ conf erenc ePisplay . py?conf Id=219381 

J. R. Espinosa, M. Muhlleitner, C. Grojean and M. Trott, "Probing for Invisible Higgs 
Decays with Global Fits," JHEP 1209 (2012) 126 |arXiv:1205.6790| [hep-ph]]. 

J. F. Gunion, Y. Jiang and S. Kraml, "Could two NMSSM Higgs bosons be present near 
125 GeV?," Phys. Rev. D 86 (2012) 071702 |arXiv: 1207. 1545] [hep-ph]]. 

P. M. Ferreira, H. E. Haber, R. Santos and J. P. Silva, "Mass-degenerate Higgs bosons at 
125 GeV in the Two-Higgs-Doublet Model," |arXiv:1211.313l| [hep-ph]. 

J. F. Gunion, Y. Jiang and S. Kraml, "Diagnosing Degenerate Higgs Bosons at 125 GeV," 
|arXiv: 1208. 18171 [hep-ph] . 

ATLAS-CONF-2012-161, "Search for the Standard Model Higgs boson in produced in 
association with a vector boson and decaying to bottom quarks with the ATLAS detector" 
(Nov 2012). 

ATLAS-CONF-201 2-160, "Search for the Standard Model Higgs boson in H -»■ t+t~ 
decays in proton-proton collisions with the ATLAS detector" (Nov 2012). 

CMS-PAS-HIG-12-039, "Search for SM Higgs in WE -»■ WWW -»■ 3£3u n (Nov 2012). 

CMS-PAS-HIG- 12-025, "Search for Higgs boson production in association with top quark 
pairs in pp collisions" (July 2012). 

CMS-PAS-HIG-12-044, "Search for the standard model Higgs boson produced in associa- 
tion with W or Z bosons, and decaying to bottom quarks for HCP 2012" (Nov 2012). 

CMS-PAS-HIG- 12-051, "Search for the standard model Higgs boson decaying to tau pairs 
produced in association with a W or Z boson" (Nov 2012). 

CMS-PAS-HIG-12-043, "Search for the standard model Higgs boson decaying to tau pairs" 
(Nov 2012). 

https : //twiki . cern. ch/twiki/bin/view/LHCPhysics/CrossSections 



22 



A. de Roeck, private communication. 

Aurelio Juste, "Standard Modek Higgs boson searches at the Tevatron " , talk at HCP2012, 
15 Nov 2012, Kyoto, Japan, 

http : //kds . kek . jp/ conf erenceDisplay . py?conf Id=9237. 

LHC Higgs Cross Section Working Group, A. David, A. Denner, M. Duehrssen, M. Grazz- 
ini, C. Grojean, G. Passarino and M. Schumacher et al, "LHC HXSWG interim recom- 



mendations to explore the coupling structure of a Higgs-like particle," [arXiv: 1209.0040 
[hep-ph] . 

M. Spira, "HIGLU: A program for the calculation of the total Higgs production cross- 
section at hadron colliders via gluon fusion including QCD corrections," |hep-ph/95 10347[ 

M. Spira, "HIGLU and HDECAY: Programs for Higgs boson production at the LHC and 
Higgs boson decay widths," Nucl. Instrum. Meth. A 389, 357 (1997) [hep-ph/9610350| . 

A. Djouadi, J. Kalinowski and M. Spira, "HDECAY: A Program for Higgs boson decays 
in the standard model and its supersymmetric extension," Comput. Phys. Commun. 108, 
56 (1998) |hep-ph/9704448] . 



F. James and M. Roos, "Minuit: A System for Function Minimization and Analysis of the 
Parameter Errors and Correlations," Comput. Phys. Commun. 10 (1975) 343. 

S. Kraml, et al, "Searches for New Physics: Les Houches Recommendations for the Pre- 
sentation of LHC Results," Eur. Phys. J. C 72 (2012) 1976 [arXiv:1203.2489| [hep-ph]]. 

F. J. Petriello, "Kaluza-Klein effects on Higgs physics in universal extra dimensions," 
JHEP 0205 (2002) 003 |hep-ph/0204067] . 

G. Belanger, A. Belyaev, M. Brown, M. Kakizaki and A. Pukhov, "Testing Minimal Univer- 
sal Extra Dimensions Using Higgs Boson Searches at the LHC," arXiv: 1207.0798 [hep-ph]. 



M. Farina, C. Grojean, F. Maltoni, E. Salvioni and A. Thamm, "Lifting degenera- 
cies in Higgs couplings using single top production in association with a Higgs boson," 



arXiv:1211.3736 [hep-ph]. 



S. Biswas, E. Gabrielli and B. Mele, "Single top and Higgs associated production as a 



probe of the Htt coupling sign at the LHC," arXiv:1211.0499 [hep-ph]. 



A. Drozd, B. Grzadkowski, J. F. Gunion and Y. Jiang, "Two-Higgs-Doublet Models and 
Enhanced Rates for a 125 GeV Higgs," |arXiv:1211.3580| [hep-ph]. 

G. C. Branco, P. M. Ferreira, L. Lavoura, M. N. Rebelo, M. Sher and J. P. Silva, 
"Theory and phenomenology of two-Higgs-doublet models," Phys. Rept. 516, 1 (2012) 
|arXiv:1106.0'034| [hep-ph]]. 



P. Teixeira-D ias, "Higgs boson searches at LEP," J. Phys. Conf. Ser. 110 (2008) 042030 
|arXiv:0804.4~l46~| [hep-ex]]. 



23 



[73] D. Choudhury, R. Islam, A. Kundu and B. Mukhopadhyaya, "Anomalous Higgs Couplings 
as a Window to New Physics," |arXiv:1212 .4652 [hep-ph]. 

[74] G. Bhattacharyya, D. Das and P. B. Pal, "Modified Higgs couplings and unitarity viola- 
tion," |arXiv: 1212.46511 [hep-ph]. 



24 



